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Abstract 

The paper deals with the study of a satellite attracted by n primary 
bodies, which form a relative equilibrium. We use orthogonal degree to 
prove global bifurcation of planar and spatial periodic solutions from the 
equilibria of the satellite. In particular, we analyze the restricted three 
body problem and the problem of a satellite attracted by the Maxwell's 
ring relative equilibrium. 

Keywords: Global bifurcation. Orthogonal degree. Restricted n-body 
problem. Ring configuration. 

1 Introduction 

The restricted n-body problem is the study of the movement of a satellite at- 
tracted by n primary bodies vifhich are rotating, at a constant angular speed, 
around an axis. Since the mass of the satellite is small, one assumes that the 
satellite does not perturb the trajectories of the primaries. We shall suppose 
that these trajectories form a relative equilibrium and, as such, are in a plane, 
let us say the {x, y)-plane. In this paper, the primaries are assumed to be point 
masses or, equivalently, homogeneous spheres. 

The purpose of this paper is to prove the existence of a global bifurcation 
of periodic solutions for the satellite, starting from the relative equilibria of the 
satellite. These solutions will form a continuum in the plane of the primaries 
and we shall also prove that there are other global branches of solutions out 
of that plane. The proof is based on the use of a topological degree for maps 
that commute with some symmetries and are orthogonal to the infinitesimal 
generators for these symmetries. We give results for a general situation and 
applications to some special cases such as the restricted three body problem 
and the Maxwell's Saturn ring, that is when there are n primaries, of the same 
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mass, forming a regular polygon, and a central larger mass, as a classical model 
for Saturn and one ring around it. However, for the general result the primaries 
may have different masses and may be located at any relative equilibrium. 

The study of relative equilibria for the restricted n-body setting is a classical 
problem and there is a vast literature for it. For instance, in the case of the 
restricted three-body problem, the local bifurcation of planar periodic orbits 
from the Lagrange points is well known (see [18], [21]). There is a huge number 
of numerical explorations for this restricted three-body problem, under a variety 
of hypotheses, such as the bifurcation near L4, where the mass of the primary 
is the bifurcation parameter above Routh's number, [3], |26] . with a period 
doubling cascade. [2B] has a study of the phase space for solutions near L4 and 
[8] treats the elliptic case where one has four periods for solutions close to L4. 
The stability of the orbits close to L4 is studied in [J and the connection from 
E3 to L4 is explored in [53]. A very complete numerical study, [5], using AUTO, 
shows the many different types of periodic orbits and the connections between 
the Lagrange points and also the secondary bifurcations along the curves in the 
x,y, fi space, where ^ is the mass of one of the primaries. From a very applied 
point of view, one may cite [15] and [11] . 

In the case of the Maxwell ring, besides the theoretical results of [17] and 
[2D], one has also many numerical results, such as [55], where the author studies 
numerically some families of solutions around the central body and around the 
ring for a low number of peripherals, with a theoretical approximation for the 
case of a satellite far from the ring. A theoretical study, with averaging tech- 
niques is given in [16] for orbits far from the set of primaries (comets) and close 
to one of the primaries (Hill solutions) . Similarly [19] proposes a regularization 
for collision orbits. Closer to the spirit of the present paper, we mention some 
of the more recent papers in the bibliography, in particular [T] , [3] , [H] and [5] , 
where a numerical classification of the different types of orbits is done. 

The paper which is closer to ours, in the sense that is based on topological 
arguments similar to ours and giving global results for the possible connections 
between the relative equilibria is [17) . for the restricted three-body problem. 

A final introductory comment, about topological methods, in particular in 
bifurcation problems, may be useful: the degree arguments, coupled with group 
representation ideas, give global information, i.e., an indication of where the 
bifurcation branches could go. Also, since the results are valid for problems 
which are deformation of the original problem, the method does not require 
high order computations and they may be applied in some degenerate cases (for 
instance it is not necessary that the bifurcation parameter crosses a critical value 
with non-zero speed; it enough that it crosses eventually). However, knowledge 
of some generic property, like a Morse condition, implies an easy application of 
the argument. This may be not the case for problems with more parameters 
(see however \T2\ ) . An immediate drawback of this approach is that topological 
methods do not provide a detailed information on the local behavior of the 
bifurcating branch, such as stability or the existence of other type of solutions, 
like KAM tori. Other methods, such as normal forms or special coordinates, 
should be used for these purposes, but they only provide local information near 
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the critical point. In a similar way, our degree arguments give only partial results 
on resonances and other tools should be used. Topological methods provide an 
interesting complement of information. 



2 Setting the problem 

Newton's equations describing the movement of a satellite, in rotating coordi- 
nates and with angular speed equal to 1, are 

x + 2Jx = W{x) with 

1 1 Tx 11^ ^ 
V{x) := ^L + ^TO,,/.«(||x-(a,-,0)||), 

where a: G is the position of the satellite and {aj , 0) is the position of a 
primary body with mass rrij . The function 4>a represents the attraction between 
the bodies, where we suppose that (j)'^ = —1/x", and we include the gravitational 
potential for a = 2. The matrices 7 and J are defined by 

I = diag(/, 0) and J = diag(J, 0), 

where J and / are the symplectic and identity 2x2 matrices. 

Here we assume that the primary bodies form a relative equilibrium. Because 
of the homogeneity of the potential, we may rescale the space so that the angular 
velocity is 1. As all relative equilibria are planar for the n-body problem, thus 
the positions of the primary bodies S must satisfy the relation 

n 

Edi — dj 
rrij- 

.•=iO¥i) ll'^^-^^ll 

The equilibria of the satellite arc just the critical points of the potential V . 
From the potential we can prove that all equilibria are planar. Now, we wish to 
find the Hessian of the potential at a planar equilibrium. 

Proposition 1 Let dj be the distance between xq = {x, y, 0) and the primary 
body (cij, 0) = {xj,yj, 0). The Hessian matrix of the potential is 

(n n \ 

I + J2mjAj,-Y,mj/d'^+'\, 

where the matrices Aj are defined by 

{x-Xj)'^ {x - Xj){y - yj) \ i_ 

' \ {x - xj){y - yj) [y-Vi? ) ^' 
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Proof. Since the function (pa{dj) has Hessian 

a + 1 [ (x-^jf {x-Xj){y-yj) 0\ j 
D^Mdj) = - Xj){y - Vj) {y - yjf - — p^-, 

\ / 

hence D'^(j)a{dj) = diag{Aj,—l/d'^~^^). From this fact we get the Hessian of V. 
■ 

Now we want to estimate the number of equiUbria provided that the potential 
is a Morse function, which is more than a reasonable condition. This is a generic 
condition, which is met in our applications, but which could not hold in some 
cases. As a matter of fact, we only need that the critical points should be 
isolated. Because all equilibria are in the plane, we may restrict the potential 
to planar points. 

Proposition 2 Let us assume that the potential of the satellite is in the plane 
with a G [1, oo). Then the potential does not have maximum points. In addition, 
if the potential is a Morse function, then 

#saddle points = n — 1 + #minimum points. 

Moreover, since the potential has a global minimum, there are at least n saddle 
points. 

Proof. The potential in the plane has Hessian D'^V{xo) = I + J2^=i ^j^j^ 
the trace of D^V{xo) is 

n 

r = 2 + (a-l)^^. (2) 

Consequently, the potential does not have maximum points as the trace is posi- 
tive for a € [1, oo). Moreover, we know that V{x) is positive and that V{x) — > oo 
as a; — > {oo, ai, a„}, then V has at least a global minimum in D,. Since the 
gradient of V{x) is dominated by the identity, for large ||x||, the critical points 
are bounded. 

Let us define the set $1 as a ball of radius p, minus small balls of radii p^^ 
with centers at aj. Since the gradient Vy points outward in dCl provided p is 
big enough, then by the Poincare-Hopf theorem the degree of WV{x) is equal 
to 1 — n. Furthermore, since V{x) is a Morse function, that is the critical points 
are non-degenerate, then this degree is the sum of the local indices. Each of 
these indices is the sign of the determinant of the Hessian matrix, that is 1 for 
a minimum and —1 for a saddle point. Then 

1 — n = degQ vy = #minimum points — #saddle points 
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3 Bifurcation theorem 



In order to explain our results, we may give a short description of the steps to 
prove the bifurcation theorem. 

We wish to remark that we follow the ideas from the book [13j , where more 
general bifurcation theorems are proven. In addition, in the thesis [5] one may 
find a systematic application to different Hamiltonian systems and situations. 

3.1 The bifurcation operator 

Our aim is to find bifurcation of periodic solutions from the equilibria of the 
satellite. First, we make the change of variables from t to t/h'. Hence, the I-k jv- 
periodic solutions of the differential equations are the 27r-periodic solutions of 

-v^x - 2vJx + W{x) = 0. 

Let be the Sobolev space of 27r-periodic functions, with the corre- 

sponding regularity. We define the collision points as the set 5* — {ai,...,a„} 
and the collision-free paths as the set 

HUW'X^) ^{x^ HU^') : x{t) ^ a,}. 

Recall that functions in this space are continuous. Hence, we define the bifur- 
cation operator / : iJ|^(M^\^') x M+ -> L^^ as 

/(x, u) = -ly^i - 2iyJx + W{x). 

In view of the definitions, the collision-free 27r-periodic solutions are zeros of the 
bifurcation operator f{x,v). Furthermore, the operator / is well defined and 
continuous. 

Now, we define the actions of the group Z2 x S*^ on i7|^(M'^\^) as 

p{k)x = Rx{t) and p{(p)x = x{t + ip), 

where R — diag(l, 1, —1) is just the reflection which fixes the plane. 

Since V{x) is invariant with respect to the reflection, the gradient W is a Z2- 
equivariant map. Moreover, since the equation is autonomous and R commutes 
with the matrix J, then 

f{p{K,ip)x) = p{K,ip)f{x). 

Therefore f{x) is a Z2 x S'^-equi variant map. 

Now, the generator of the group on the space is D^{p{'p)x)^p^Q = x. 
As the operator f{x) satisfles the equality 

(/W, = (-^' |i|V2 + = 0, 

then the operator f{x) is orthogonal to the generator x in L"^^. Given this 
condition we say that the operator f{x) is a Z2 x S'^-orthogonal map. The 
orthogonality corresponds to the conservation of energy. 

Finally, since all the equilibria are planar, the isotropy subgroup of an equi- 
librium xq is Z2 X S"^. This means that all equilibria are fixed by the action of 
Z2 X S^. 
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3.2 The Lyapunov-Schmidt reduction 

Wc want to use the orthogonal degree in order to prove bifurcation, but since 
this degree is defined only in finite dimensions, we need to reduce the bifurcation 
operator to finite dimensions. To achieve this, let us set the Fourier series of 
the bifurcation operator as 

fix) = J2 {l^'^^^i - '^il'^J^i + 9i) e'", 

where xi and gi are the Fourier modes of x and \/V{x), respectively. Since 
Z^i/^J — 2«/i^J is invertible for all Vs except a finite number, we can make a 
Lyapunov-Schmidt reduction to a finite space. In fact, we perform a global 
reduction, using the global implicit function theorem, with the right bounds 
taking care of the collision points 

In that way, we get that the zeros of the bifurcation operator are the zeros 
of the bifurcation function 

f{xi,X2{xi,v), = ^ [fv^xi - 2ilvJxi + 
\i\<p 

where xi corresponds to the 2p + 1 modes and X2 to the complement. 

Consequently, the linearized bifurcation function at an equilibrium xq is 

f\xo,v)x-^ = {I'^v'^I -'i.ilvJ + D'^V{xq)) xief'^* . 
\i\<p 

In fact, '^V{x) = D'^V{xq){x — x^) + .... close to .tq and the Fourier compo- 
nents oi X — xq are xi for Z ^ and we rename the stationary mode as xq. 

So the linearized bifurcation equation has blocks M{lu) for I e {0, 
where the block M(A) is 

M(A) = - 2iXJ + D'^V{xo). 

3.3 Irreducible representations 

In the following part, we analyze the symmetries of the group Z2 x S^. Since 
the action of {k, G Z2 x on Fourier modes e''*a;; is 

piK,^){e'''xi)=Re''^{e''*xi), 

then the action on the block M(li') is given by p{k, ip)xi — _Re*''^a;(. 

Now, as the action of Z2 on is p{k) — diag{l,l, —1), the space has 
two irreducible representations: Vq = x {0} and Vi = {0} x C. That is, the 
group Z2 acts on Vo as p{k) — 1 and on Vi as p{k) = —1. Hence, by Schur's 
lemma we know that the matrix M(A) must satisfy 

M(A) =diag(Mo(A),Mi(A)). 
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Actually, from the explicit Hessian D'^V{xq) we have 

n 

Mi(A) = - ^ mj/d'^+^ and (3) 

Afo(A) = X^I - 2iJX + (^I + Yl "^J^J" j ■ 

Consequently, the action of the group Z2 x on the block Mo{i/) is (k, 1^9)2; = 
e^'^'x. Therefore the element (k, 0) leaves fixed the points for Mo(z^), so the 
isotropy subgroup for MQ{iy) is the one generated by (k, 0), 

Z2 = ((k,0)). 

For Afi(:^) the action of the group Z2 x 5^ is {K,ip)x — —e^'^x. It follows that 
(kjTt) leaves fixed the points for Mi{i'), thus the isotropy subgroup for Afi(j^) 
is generated by (k, tt), 

= ((At,7r)). 



3.4 The orthogonal degree 

The orthogonal degree is defined for orthogonal maps that are non-zero on the 
boundary of some open bounded invariant set. The degree is made of integers, 
one for each orbit type, and it has all the properties of the usual Brouwer degree. 
Hence, if one of the integers is non-zero, then the map has a zero corresponding 
to the orbit type of that integer. In addition, the degree is invariant under 
orthogonal deformations that are non-zero on the boundary. The degree has 
other properties such as sum, products and suspensions, for instance, the degree 
of two pieces of the set is the sum of the degrees. 

Now, if one has an isolated orbit, then its linearization at one point of the 
orbit xq has a block diagonal structure, due to Schur's lemma, where the isotropy 
subgroup of xq acts as Z„ or as S^. Therefore, the orthogonal index of the orbit 
is given by the signs of the determinants of the submatrices where the action 
is as Z„, for n — 1 and n — 2, and the Morse indices of the submatrices 
where the action is as S^. In particular, for problems with a parameter, if 
the orthogonal index changes at some value of the parameter, one will have 
bifurcation of solutions with the corresponding orbit type. Here, the parameter 
is the frequency v. 

Any Fourier mode will give rise to an orbit type (modes which are multiples 
of it), hence one has an element of the orthogonal degree for each mode. Fur- 
thermore, if x{t) is a periodic solution, with frequency v, then y{t) = x(nt) is 
a 27r/n-periodic solution, with frequency v/n. Hence, any branch arising from 
the fundamental mode will be reproduced in the harmonic branch. If one wishes 
to study period-doubling, then one has to consider the branch corresponding to 
TT-periodic solutions. 

The complete study of the orthogonal degree theory is given in jl3) . 
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Theorem 3 Supposing that the matrix M(0) — D'^V{xq) is invertible, we de- 
fine 

r]k{X) = <7{nk{\- p) -nk{\ + p)), (4) 

where a — sgn(det Afo(O)) and 7ifc(A) is the Morse index of Mk{X) for k £ {0, 1}. 

In general, if xq is an isolated critical point, then a is the index of W{x) 
at xq. 

Ifr]k{h'k) is nonzero, then the equilibrium has a global bifurcation of periodic 
solutions starting from the period 2TT/i'k with isotropy group Gk- 

Proof. Since Mi(0) is a negative number, the sign of the determinant of M(0) 
is the opposite of a. Furthermore, there will be a change of the Morse number 
only at values of A where Mi (A) is or where the self-adjoint matrix Mo(A) has 
one of its two eigenvalues equal to ( the other is not 0, given that the trace 
is positive). Finally, since A = liy, what happens for the fundamental mode 
(/ = 1) is reproduced for higher modes and frequencies which are quotients 
of the fundamental frequency by the mode I. Here we take the fundamental 
mode. One is then in the position of applying Proposition 3.1, p. 255 of [T5] . 
after one sees the change of orthogonal index. Finally, if xq is an isolated 
critical point, then one may perform an orthogonal deformation of the map to 
{\/V{x), M{li')xi), for I € {!,..., p}, near (xo,!'*;), with a jump at Vk given by 
the above formula. ■ 

We say that the bifurcation is non- admissible when either: i) the global 
branch goes to infinity in norm or period or ii) the branch ends in a collision path. 
In any other case we say that the bifurcation is admissible. By global bifurcation 
we mean either that the bifurcation is non-admissible or, if the bifurcation is 
admissible, that the bifurcation branch returns to other bifurcation points and 
that the sum of the jumps of the indices at the bifurcation points, rjkii'k), is 
zero. 

4 Spectral analysis 

Now, we wish to find the bifurcation points of an equilibrium. In order to do 
so, we need to analyze the spectrum of the blocks Mo(A) and Mi (A). But let 
us first find the symmetries of the solutions that bifurcate from these blocks. 

For Mo (A) we get solutions with isotropy subgroup Z2. As k S Z2 has action 
p{K)xa{t) = Rxo{t), this means that the solutions with symmetry Z2 satisfies 
xo{t) — Rxo{t), i.e. z{t) — 0. Therefore, solutions with symmetry Z2 are just 
planar solutions. 

For Ml (A) we get solutions with isotropy subgroup Z2. As {k,tt) G Z2 has 
action p{{K,Tr)xo{t) — RxQ{t + tt), then the solutions with symmetry Z2 satisfy 
Xo{t) — Rx^it + tt), i.e. 

x{t) = x{t + tt), y{t) = y{t + tt) and z{t) = -z{t + tt). (5) 

Since the projection of this solution on the (x, j/)-plane is a 7r-periodic curve, that 
solution follows twice this planar curve, one time with the spatial coordinate 
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z{t) and a second time with —z{t). Consequently, there is at least one to where 
z{to) = z{to + n) =0. For instance, if only one to exists, then the solution 
looks like a spatial eight near the equilibrium. For this reason, we will call 
eight-solutions the solutions with isotropy subgroup Z2. 

Remark 4 Actually, the solutions of the satellite are defined in rotating co- 
ordinates, so that the periodic solutions are in general quasiperiodic in fixed 
coordinates. 



4.1 Planar solutions 

Let T and D be the trace and determinant of the matrix Afo(O). We point 
out that the block Mo(0) is just the Hessian of the planar potential at the 
equilibrium point. In addition, in the first section we have proven that the 
trace T is always positive. Now, we want to show that the bifurcation depends 
essentially on the sign of D. 

Proposition 5 Let us define u± as 

1/2 



u± = (2- T/2 ± ^7(2 - r/2)2 - 



(a) If D <0, then xq has a global bifurcation of periodic planar solutions from 

2n/i'+ with 

rjo{j^+) = -1. 

(b) If D > 0, (2 - T/2)2 > D and T < 4, then xq has a global bifurcation of 

periodic planar solutions from 2^/1/+ and 2Tr/u- with 

T]o{i'+) = 1 and rjoiv-) = -1. 



Proof. Since Mo(0) is selfadjoint, there is an orthonormal matrix P € SO{2) 
such that A'/o(0) = P^AP, where A is the eigenvalue matrix diag(Ai, A2). Since 
Mo(z/) = J — 2i Jz/ + Afo(O) and J commutes with P, then 

PMo{u)P^ = diag(z/2 + Ai, z/2 + A2) - 2u{iJ). 

In view of T = Ai + A2 and D = X1X2, the determinant of Mo{u) is 

det Mo{u) = v^- 2(2 - T/2)v'^ + D. 

It follows that the determinant has the factorization 

detMo(i/) = (i/2 - vl){v^ - 1/2). 

Consequently, the Morse index of Mo{v) can change only at ±,Jv±. 

For (a), only is positive, and a — sgn(_D) — —1. The Morse index of 
Mo(0) is no(0) = 1 due to D < 0, and no(oo) = due to the fact that Mo{v) has 
only positive eigenvalues for u big enough. Therefore r?o('^+) = cr(l — 0) = — 1. 
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For (b), both numbers v± are positive, and a = sgn{D) = 1. Moreover, 
we see that the determinant of Mq(v) is negative between z/_ and z^+, thus 
no(i^) = 1 for i> G As the Morse index at infinity is 710(00) = 0, we 

conclude that 77o(i^+) = 1 — 0. Now, the Morse index of Mo(0) is no(0) = 2 if 
T < and no(0) = if T > 0. It fohows that r]o{iy^) = 2 - 1 if T < and 
,yo(,y_)=0-lifT>0. 

Note that this proof is independent of the form of the potentiaL For the 
case of the specific potential of this paper, equation (2) implies that T > 0. ■ 

Remark 6 In the case (b), the two local bifurcations can locally collide when 
the resonance condition — mv^ holds. Moreover, it is easy to prove that the 
resonance condition is equivalent to 

(4 - T)D-^/^ e {to + m-^ : m e N}. 

Remark 7 In all other cases different from (a),( b) and (2 — T/2)^ — D, 
there is no bifurcation, if D is not 0, since then the matrix Mq{\) is always 
invertible. In addition, in the case (2 — T/2)^ — D > 0, both points v± coincide 
and rjQ{v^) = 0, then we cannot assure or discard the existence of bifurcation, 
but probably of a different kind, as found in ^ and \26f . Finally, if D = 0, then 
= and = (4 — TY^^ , if T < 4, i.e. V{x) is not a Morse function at xq. 
In this last case, one may have a bifurcation of relative equilibria if the masses 
of the primaries are chosen as parameter and one has a change in a, when one 
of masses crosses the critical value, or one could have a secondary bifurcation 
of periodic solutions if the unfolding has the right properties, see \12f . However, 
in the applications of the present paper, the potential is a Morse function. 

Remark 8 Actually, the satellite equation on the plane is a Hamiltonian system 
with two degrees of freedom. We can relate the linear stability of the system 
with the bifurcation analysis. Indeed, it can be proven that the equilibrium xo 
is linearly stable on the plane if and only if condition (b ) is satisfied. Note that 
one could argue about the usefulness of a bifurcation result for the satellite if the 
arrangement of the primaries is unstable. This is a quite valid argument from the 
practical point of view, taking into account the reality of this model for a problem 
of mechanics. However, the mathematical result is independent of the stability of 
the primaries and furthermore, as proved in 0/ and in an article in preparation, 
the primaries may loose their stability and generate stable periodic solutions of 
the whole system. In that case, it is much simpler to prove the bifurcation of 
periodic solutions for the satellite, assuming, as a first approximation, that the 
primaries are at their position of relative equilibrium. Hence, the mathematical 
study of the bifurcation is also justified in this framework. 

In the case of the Maxwell ring, it is well known that the system of the 
primaries is unstable if n is between 3 and 6 and the stability is treated, for n > 6 
and large central mass, in \24^ , 128}/ and others. A complete mathematical study 
of the stability is given in flOf . Thus, if one insists, on physical grounds, that 
the stability of the relative equilibrium configuration must be insured in other to 
have a study of the bifurcation, one has to restrict to the case n > 6 and large 
central mass, or assume that the primaries are fixed in the rotating frame. 
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Remark 9 Because there is only one bifurcation value for the frequency in case 
(a), the global branch cannot return to the same equilibrium point, so the bifur- 
cation branch is non-admissible or it is connected to the bifurcation point of 
another equilibrium. In fact, if the potential is a Morse function, then one 
should get a connection to the small period branch of a minimum ( that is with a 
jump of 1 ). This implies that, in this case, there are at least n—1 non admissi- 
ble branches starting from saddle points, (see our previous proposition). In JJ7| /, 
one finds other possibilities for branches starting from a minimum, and v^, for 
the restricted three-body problem. 

4.2 Spatial solutions 

As before, let v± be the points where Mo(A) is not invertible . 
Proposition 10 Let us define vi as the positive root of 

n 

Then every equilibrium xq has a global bifurcation of periodic eight solutions 
with 

r]i{vi) = a. 

In addition, the local bifurcation branch from 2t:/ui is truly spatial, z{t) ^ 0, 
provided the nonresonant condition vi ^ v±/2l holds. 

Proof. It is clear that Mi{v) is zero only for ±z/i. Since Mi(oo) is positive and 
Mi(0) is negative, the Morse indices at infinity and zero are ni(oo) = and 
ni(0) = 1. Therefore = (t(1 — 0). Thus, one has the global bifurcation of 

periodic eight solutions. 

It remains only to prove that the solutions are truly spatial. In order to 
achieve this, we need to prove the nonexistence of solutions of the kind 

x{t) = x{t + tt), y{t) = y{t + vr) and z(t) = (6) 

near (xq, vi). In fact, the solutions are in the fixed point space of the group 
Z2 X Z2 generated by k G Z2 and tt ^ . 

Now, the restriction of the derivative of the bifurcation equation to the fixed 
point space of Z2 x Z2 has blocks Mq{21vi). Since the matrix Mq{v) is invertible 
except for the points v±, and we suppose v± ^ 2lvi, the blocks Afo(2/i'o) are 
invertible. Consequently, the derivative of the bifurcation equation in the fixed 
point space of Z2 x Z2 is invertible. Therefore, we get the nonexistence of planar 
solutions ^ near {xq^vi) from the implicit function theorem. ■ 

Remark 11 Although the nonresonant condition vi 7^ v±/2l is sufficient to 
assure that the bifurcation from 27: /vi is really spatial, it is not a necessary con- 
dition. If one considers the full three-dimensional problem, without any special 
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symmetry (except periodicity), then, if one has the resonance v± — 2lvi, the 
jump of orthogonal index has two components 771(1^1) for the fundamental mode 
and ?7o(j^±) for the 2l-mode. Since this jump is different from the one caused 
by the rescaling of the jump for the solution in the fixed-point subspace 0/Z2, 
which has only the second component for the 2l-mode, one obtains a new branch 
of periodic solutions. In the case of the restricted three-body problem, this is the 
branch given in JIT] /. 

5 Applications 

5.1 A Morse Potential 

We have proven that the potential for the satelhte problem has at least n saddle 
points and a global minimum, provided it is a Morse function. Consequently, 
we get the following result: 

Theorem 12 Each one of the saddle points has a global bifurcation of planar 
periodic solutions and a global bifurcation of periodic eight solutions. 

Each of the minimum points satisfy one of the following options: (a) it has 
two global bifurcations of planar periodic solutions and one bifurcation of periodic 
eight solutions, or (b ) it has one bifurcation of spatial periodic eight solutions. 

For the planar bifurcation, each saddle point has a bifurcation with index 
770 = —1 and each minimum point has two bifurcations, if any, one with rjQ = 1 
and another with 770 = —1- Because an admissible bifurcation branch has sum of 
indices ryo equal to zero, the sum over all admissible branches is 0. If Sa denotes 
the number of saddle points which belong to an admissible branch, m_a the 
number of minima with jump of —1 which are on an admissible branch and 
m+a those with jump 1, one has that Sa + m-a = rn+a- Let Si, m^i, m+i be 
the numbers of points which are on non-admissble branches and let s be the total 
number of saddle points, m the number of minima (including mg those which are 
not on any branch), then one gets that m — mQ-^m^a + rn^i = mo +TO+a + 7n+i 
and, since s = n — 1 + to, one has Si + m^i — m+i = ri — 1 + m, that is the 
number Si + m-i of points with jump —1 belonging to non-admissible branches 
is at least n — 1 + ?Ti. Thus, the number of points on non-admissible branches is 
at least the number of saddle points. 

Now, since every minimum has a spatial bifurcation with 771 = 1 and every 
saddle point has a spatial bifurcation with rji = —1, then a bifurcation branch 
of eight solutions is non-admissible or the total number of saddle and minimum 
points that it connects is the same and the number of saddle points which are 
on non- admissible branches of eight solutions is at least n — 1. 

5.2 The restricted three-body problem 

In the restricted three-body problem, the primary bodies are at ai = (1 — /i, 0) 
and 02 — {—fi,0) with masses toi = fi and 7712 = 1 — fi. Hence, the potential of 
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the satellite is 

1 - ^ 
nx) = -||/xf +^m,</>4lk-(a„0)||). 

j=i 

This problem is well known on the plane, see for instance [21] • There are 
only five equilibrium points called Lagrangians. Two of these equilibrium points 
form an equilateral triangle with the primary bodies ai and 02, and they are 
minima of the planar potential. Three of the equilibrium points are coUinear 
with the primaries, also called Eulerian points, and they are saddle points of 
the potential. All of these relative equilibria are non-degenerate, that is V{x) 
is a Morse function. 

Also, it is well known that the minimum points have two bifurcation fre- 
quencies for fi < /ii, where fii = (1 — (a + l)^^ a{30 — a) — 33)/12)/2, 
when a is in the interval(15 — 8-\/3, 3), is the critical Routh ratio and without 
any restriction on if a belongs to the interval (1, 15 — 8-\/3). This comes from 
the fact that the trace T — a + 1 and the determinant D — 3{a + 1)^m(1 ^ /4, 
with the conditions T < 4 and {2 - T/2)'^ > D. In that case, 

1^1 = (3- a± v/(3 - a)2 - 3(a + l)2/i(l - /2. 

Note that vj^jv^ tends to infinity when [i, tends to 0, thus there is an infinite 
number of resonance values for /x, when /i goes to 0. 

For the saddle points we have only the bifurcation point where 

v\ = \-{a^ \)vll2 -f ((a + 1)2^.4/4 _ 2(a - \)vlYl^ 

with v\ — ^'^j/'^j > Ij since D'^V{xq) — diag{l + aiyf, 1 — i^f). 

Consequently, we get the classical global bifurcation of planar periodic solu- 
tions, with at least three equilibria on non-admissible branches, see |17j for the 
case a = 2. 

Now, we wish to find bifurcation of spatial periodic eight-solutions. 

Theorem 13 In the restricted three-body problem each one of the five equilibria 
has a global bifurcation of spatial periodic eight-solutions. 

Proof. We only need to prove the nonresonant condition 1^1 > v±/2l at equi- 
librium points. For the triangular Lagrangian points we have that vi — 1 and 
v± £ (0, y/{3 — a) for e (0, 1), therefore i/i > v±/2l. 

For the collinear Lagrangian points, since vj^ is given in terms of vi, we need 
to prove that Af'vl ^ v\, or equivalently (a + lYvt/^ - 2{a - l)vl ^ ((4^^ -f- 
(a — l)/2)i'^ — 1)2. The last inequality is also equivalent to avf — 2bvi +1^0, 
where a = (4P + (a - l)/2)^ - {a + 1)2/4 and b ^ 4P - {a - l)/2. But since 
b^-a = {a + 1)2/4 - 8(a - 1)P < is satisfied for ah 1, if a £ (15 - 8^3, 3), 
then the quadratic equation aiyf — 2bvl -1-1 = does not have solutions and 
APvI ^ v\. On the other hand, if a e (1,15 - 8^3) and I = 1, then the 
quadratic equation has its largest root less than 1, which contradicts the fact 
that at the saddle point vi > 1. Thus, there is no resonance and the branch is 
truly spatial and at least one branch is non- admissible. ■ 
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5.3 The Maxwell's Saturn ring 



In this section, wc analyze the satelhte problem when the primaries form a 
polygonal relative equilibrium. Hereafter, we identify the real and complex 
planes. 

The polygon consists of one body of mass /i at oq = 0, and n bodies of 
mass 1 at each vertex of a regular polygon, for instance aj = ae*-''- for j G 
{l,...,n}, where = 2iT/n. It is easy to prove that the positions aj form a 
relative equilibrium provided that a"+^ = s + fj,, where s is defined by 



2«^sin«-i(jC/2)' 

Moreover, we can make the change of variable x = au in such a way that 
the equation is u + 2Ju = W{u) with the potential 

V{u) = - ||7«f + -—M\u- [e'^^M) + -^M\\u\\). 

Now we point out that the case n = 2 with /i = is just a particular case of 
the restricted three-body problem, hence we shall analyze only the cases n = 2 
with > and n > 3 with /U > 0. 



Existence of equilibria 



Remember that all equilibrium points of the satellite are in the plane. So, we 
assume, for this purpose, that the satellite is in the plane, i.e. the potential is 



1 1 



!)■ 



s + n 



with u e 



Proposition 14 For i-i = 0, the origin uq = is a critical point. In addition, 
we have forn>3 that D'^V{0) = XI with A > 0. 

Proof. That the origin is a critical point follows from the fact that 

n 

V„y(0) = -ye^^«=0. 
s 

Now, since D'^V{0) has real eigenvalues and D^V{0) is £)„-equivariant, by 
Schur's lemma we have D^V(O) = AI for n > 3. That A > is due to the 
fact that the trace T = 2A is always positive. ■ 

Now for u 7^ 0, we may simplify the analysis if we change to polar coordi- 
nates. For these coordinates the potential is 



V{r, ^) = rV2 + 



n 



j{j(-<p) 
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Let us observe that the potential V is D„-invariant for the action p{C)u = 
e^^u and p{k)u = u, thus, critical points will be D„-orbits of points. It follows 
that the potential V{r, (p) is even and 27r/n-periodic in (p, hence, we may restrict 
our analysis to points with ip G [0, Tr/n]. 

Now, we will show that the potential has three orbits of critical points. To 
achieve this goal, we need first to prove the following lemma. 

Lemma 15 For n > 3, the derivative Vr at e"/" is negative, 

Vril, n/n) < 0. 

Proof. The derivative of V{r, (p) is 



Therefore, at e"/", we have 



+ /i \^^4^2"sin"-i(j- 



+ y—'i 2" sin""^ (j - l/2)C/2 ; s + fi' 

where a is the sum between parentheses. 

So it remains to prove that s < a. In order to do so, we need some inequal- 
ities. Since n > 3, we have the first inequality 

n-l 

2"5 = V - < 2 V = 

where equality holds for n odd. Similarly, we have the second inequality 

1 



2v-y — ^ >2 y 

^ sin"-i(7 - l/2)C/2 - . 



where equality holds for n even. Finally, since sm{j — l/2)C/2 < sinjC/2 for 
j G [1, n/2], then we have the third inequality 



> 



sin"-i(j-l/2)C/2 sin"-i(jC/2)' 



The fact a > s follows from these inequalities. ■ 

In [5], one may find an integral representation which is used to prove the 
next corollary. In addition, a direct proof of the integral representation and of 
this corollary will be given in the last section. 

Corollary 16 For a £ (1,3), the derivative Vr{r,p) is the product of — sininip) 
with a positive function uj{r, p), 

Kp(r, ip) ^ - fi\n{nip)uj{r, ip). 
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We may now prove the existence of Z„-orbits of equilibrium points. 

Proposition 17 For a € (1,3) and ri > 3 there are three orbits of critical 
points. We are showing only the points of the Zn-orbits with ip € [0,7r/n]; 

(a) If fJ- (z (0, oo), there are two saddle points at r2 and ri, with < r2 < 1 < ri, 

and there is a minimum, point at rse''^/", with r-^ > 1. 

(b) If ^ = 0, there are two saddle points at ri and r2e^^^^, with < r2 < 1 < ri, 

and there is a minimum point at r-^e}'^!'"- , with > 1. 

Furthermore, there are no other critical points when (p € [0,7r/n). 

Proof. Since V^{r,(fi) = — sin{n(p)uj{r,(p), with a positive function co{r,(p), 
then V^{r,ip) = only for ip = kn/n. Furthermore, at these points we have 
V^^{r,kn/n) = —noj{r,ip)cosk'jT. Consequently, the critical points must be in 
(p e {0, 7r/n} with 

V^^{r,0) < and V^^{r,n/n) > 0. 

Thus, in order to find critical points, we need to look only for points where 
Vr{r,ip) = 0, with (p = 0,7r/n. 

Before we start finding critical points, we wish to prove that all the critical 
points with ip = are saddle points. The trace of D'^V{xo) at a critical point is 

Similarly, it is easy to see that the determinant of D'^V{xq) at a critical point is 

D = VrrV^^r-^. 

Now, since T is always positive and V^^(r,Q) is always negative, then Vrr{r,0) 
is positive. Consequently, all critical points, with p = 0, satisfy 

Vrr{r,0) > and V^^{r,0) < 0. 

For ji G [0,oo), the potential V{r,0) goes to infinity when r {l,oo}. 
Hence, the potential has a saddle point at ri G (l,oo). Now, if there were 
another critical point r* in (l,oo), then V^r(?"*,0) would be positive. In that 
case there would be another critical point between ri and with Vrr{r, 0) < 0. 
But that cannot happen, and consequently ri is unique in (l,cx)). 

For /i e (0, oo), the potential V{r, 0) goes to infinity when r — >■ {1, 0}. Hence 
the potential has a saddle point with r2 € (0)1)- As before with ri, we can 
prove that r2 is unique in (0, 1). 

For 1-1 = 0, remember that K-(0, = and Vrr{0, ^) > 0. Then, by a similar 
argument to the uniqueness of ri we can prove that the potential V{r, 0) does 
not have critical points in (0, 1). Now, for (p — 7r/n, since V^(l, 7r/n) is negative 
and K-(0,7r/n) = with T4r(0,7r/n) > 0, there must be a r2 < 1 such that 
Vr(r2,7r/n) = with lrr(^'2,7r/n) < 0. Consequently r2e*'^/"' is a saddle point. 
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For fi S [0,00), since 1^(1, 7r/n) is negative and since V^(r, 7r/n) goes to 
infinity as r — > 00, there is a critical point € (1, 00) such that Vrrir^, ix jn) > 0. 
Therefore rse*'^/" is a minimum. ■ 

In the article [3] , the existence of these three orbits of equilibrium points is 
proven, as well as their stability. However, our proofs are simpler. 

For n = 2 and > we can prove the previous proposition with the same 
argument, except for the existence of rj,. Instead, we get the existence of a 
ra G (0,00) because the potential ^(r, 7r/2) goes to infinity when r — >■ 0, 00. 

Now, in [3], the question of the existence of more critical points was left 
open. Actually, for n = 1 and > we can prove the following: 

Proposition 18 For n — 2 and fi > the previous proposition is true and 
there are no other critical points. 

Proof. It remains only to prove that r^ is in (l,oo) and is unique. Let us 
define /(r) = — 2(r^ + After some computations, we find that the 

derivative Vr{r, 7r/2) satisfies the equality 

(s + ^l)Vr - r(/(r) + s) + ^i{r - r""). (8) 

Let us denote the /^-dependence of the potential as V{r, /i). Therefore, from 
the equality dH), we have that Vr{r, 7r/2; /i) < Vr{r, 7r/2; 0) for r < 1. Now, as the 
three body problem is the case n = 2 with /i = 0, we know that Vr{r, 7r/2; 0) = 
only at the triangular Lagrangian point r = V2- Furthermore, V^r('', 7r/2; 0) < 
for r < 1, and hence Vr{r, 7r/2; /i) < for r < 1. 

Now, let us analyze the case r > 1. From ([S]), we see, for the second deriva- 
tive, that 

(5 -f y)Vrr = [rf + f ) + s + ^i{l + ar-("+i)). 

Since rf'+f — 2 (r^a — l) [r^ + l) '■"^^^^^ jg ^ positive function, then Vrr{r, 7r/2) > 
for r > 1. From this statement, we conclude that y,.(r, 7r/2) has only the crit- 
ical point ra e (l,oo). ■ 

We proved that there may be more critical points only ii ip — n/n. And 
indeed, for n > 3 we can find more critical points when /i is near zero. 

Proposition 19 For n > 3 and fi near zero the potential has also a minimum 
and a saddle point at r4e''^/" and r^e^^^^ with r4 < r^ < 1. On the other hand, 
for ji large, rge^'^/" is the only critical point on that line. 

Proof. As before, we represent the dependence of the potential in /i as Vr{r, ii/n] /i). 
Remember that K-(0, (/3;0) = with Kr-(0,V';0) = A > for n > 3, then there 
is a r* e (0,£) such that Vrir^,,!! /n;Q) > 0. Therefore, Vr{r^,Tr/n; fi) > for 
fi near zero due to the continuity. Gathering data, we get Vr{0,T:/n) — ~oo, 
Vr{r,,,Tr/n) > and VriljTr/n) < for near zero. Consequently, there are 
two points r4 and rs where Vr{r,Tr/n) is zero with r^ < r^ < 1. Moreover, the 
second derivative satisfies Vrr{r,TT/n) > for r close to r4 and Vrr{r,Tr/n) < 
for r close to rr,. Therefore, r4e'^/" is a minimum and r^e^^^^ is a saddle point. 
On the other hand, for n large it is easy to see that Vr is strictly increasing. ■ 
The existence of the solutions r4e"/" and rse"/" was pointed out in the 
paper [I]. 
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Existence of bifurcation 



At the saddle points we have the fohowing result: 

Theorem 20 The potential has two "Ln- orbits of saddle points for n > 2, and 
one more when n > 3 and fi is near zero. Furthermore, each one of the sad- 
dle points has one global bifurcation of planar periodic and one bifurcation of 
periodic eight solutions. 

Proof. The saddle point on the line iy9 = is non-degenerate, while the critical 
points on the line if — ir/n are isolated, since Vr is locally analytic. Hence the 
index at r^ will be —1, unless r^ and r^ coincide, in which case the index would 
be 0. ■ 

Also at the orbit of minimum points we have the following: 

Theorem 21 The potential has one Tin-orbit of minimum points for n > 2, and 
one more when n > 3 and n is near zero. Moreover, provided fi is big enough, 
each minimum point has two global bifurcations of planar periodic solutions and 
one global bifurcation of periodic eight solutions. On the other hand, if a>2 and 
/i is small, the minimum v/^e^'^l'^ has no bifurcation of planar periodic solutions 
and it has a global bifurcation of spatial eight solutions. 

Proof. Since the minima are isolated, with a — I, we only need to confirm 
that the bifurcation condition (b), T < 4 and (2 - T/2)'^ > D > 0, is satisfied 
at Tac*'^/" provided that fi is big enough. 

As ra is a critical point, i.e. Vr{r^e'^^ fi) = 0, from ([7]) we can see that 
r3(^) — )■ 1 when /x — )■ ck). From the definition ([ij of Aj, the matrix 



Given that T{^i) a + l and D^p) ^ when ^ oo, then (2 - T/2f - D ^ 
e > for a £ (1,3). Consequently, for a € (1,3), at the minimum point the 
bifurcation condition (b) holds provided fi is big enough. Finally, for the minima 
inside the unit disc, one has that di and d„ are less than 1, hence, for a > 2 
one has that T > 4. ■ 

Remark 22 As a consequence of the previous proposition, we get that the min- 
imum point r^e^^/^ is linearly stable for fj, big enough. This is one of the aims 
of the article J&j where the stability, for the system of the primaries and the 
satellite, is proven for n > 7 and fi big enough. 




converges, when fi oo, to the matrix 




(cos 7r/n)^ 
C0S7r/nsin7r/n 



cosTr/nsinvr/n 
(sin7r/ri)^ 
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Remark 23 For 71 > 3 with fj. = 0, as we have seen before, at the origin xq = 0, 
we have Mo(0) — XI. Actually, since the trace T > 2, we can prove that the 
condition for bifurcation (b) is not satisfied. Hence the origin is a minimum 
point without bifurcation of planar periodic solutions. 

On the other hand, the origin is a minimum with one bifurcation of spatial 
eight periodic solutions. Moreover, we can prove, from the symmetries, that the 
bifurcating solutions satisfy x(t) — 0, y{t) = and z{t) ~ —z(t + tt). In fact, 
we can find z{t) by quadrature from the equation z ~ W^z), with V{z) := 
■^(paW -2^ + 1); with close to n. Recall that, at it is well known, that in this 
case the system of the primaries is linearly unstable. 

Remark 24 The study of the bifurcation of periodic solutions, in the plane and 
also in space, for the full system of primaries, will be published in another paper. 



6 Appendix: Integral representation 

Let us define the sum S{r, ip) as 

n 



In [3] the following integral representation of S{r, (p) is proved. We shall give 
here a direct proof using Cauchy integrals. 

Lemma 25 For (3 e (0, 2) and r e (0, 1) the function S{r, ip) has the integral 
representation 

S{r,p) ^ -s\n{TTP/2) j j— ^—jj^ fir) dr with 
fir) = I 



2 ■ 



We are now in a position of proving the corollary on Vip. 
From the integral representation we get that 



. . - . r 1 2(rT)" 1 - {rrf 

= — sui(no9) — sm— — / riTTT;—, „ , ,„ jdr 

^ ^ - 2 io {T-^-\fl^T{\-r^TYI^\\\^(Tre-'v)nf 



for r e [0, 1). Therefore S^(r, p) is the product of — ^mirvp) with the function 
between parentheses, which is positive. Moreover, since the sum S(r) satis- 
fies the equality 5'(l/r) = r^S{r), we conclude that Sip{r,ip) is the product of 
— sin(nip) with a positive function for (3 G (0, 2) and r ^ 1. 

We use this integral representation to prove the following: Set /S — a ~ 1, 
then we have (pair) = l/{(3r^). Now, we can express the potential V{r,ip) in 
terms of S{r, ip) as 

V{r, ^) = rV2 + ^Mr) + -^ks{r, p,). 

s + /i s + /i p 
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Since V depends on (p only through S(r,(p), wc conchidc that V^{r,(p) is the 
product of — sin(n<^) with a positive function for a = /3 + 1 G (1, 3). 

We may now prove the last lemma: 
Proof. Let us define the function 'w{z) as 

1 

•w{z) 



[^-l-l]/5/2' 



This function has an analytic extension to C— [0, 1]. Indeed, using the principal 
branch of the logarithm we can extend it as 

y,^^,) = g-(/3/2)[log|z-i-l|+iarg(z-i-l)l_ 

Let w^(r) be the limits w^(r) = lim£_>.o 'w{r ±i\e\) for r G (0, 1), then 
u;+(r) = e"'^"' _ „ , and w~{r) = 



(r-l _ l)/3/2 " " ^ V' ; _ ^^;3/2- 

Let be the set of points 

n^ = {\z\<l/s: |^-r| >£forr e [0,1]}. 

As the function w{z)f{z) is of order 0(1/2^+^/^) when z oo, if /3 > then 
the integral over the circle of radius 1/e goes to zero when e — )■ 0. Moreover, 
since the product w{z)f{z) is of order 0{z^/'^~^) when 2 — >• and of order 
0((1 — z)~^/'^) when z ^ 1, then for /3 < 2 the integrals over the half circles 
around = and = 1 go to zero when £ — )• 0. Consequently, we have that 

lim / w{z) f {z)dz = / [w^ {t) — w~ {T)\f {T)dT 
Jdn, Jo 

= (e-'/3--l) [\-i^r)f{T)dT. 
Jo 

Now, the function ■w(z)f{z) has n poles in C — [0, 1] and another one at 
z = r^2, but the residue at z = r"^ is zero because (3/2 e (0, 1). The other n 
poles are the roots of the polynomial function 

giz) = ||l - (zre-*^)"f = 1 + (rz)^" - 2(rz)" cosnv?. 

Consequently, the poles are found at the points z~^ — (re~*'^)e*-'^ for j = 
0, n — 1. As (TZj)"^ — e™'^, the derivative of g at the pole zj is 

g'{zj) = 2/1^7 ie*"'^(e"'^ - cosn<^) = 2inzJ^e'"'^ sinmp. 

Consequently, the residue of w{z)f{z) at the pole Zj is 

1 1 - e^™"^ 



res^w{z)f{z) 



- 1)(1 - rHj)Y/^ z^g'iz^) 
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Moreover, since r'^zj = Zj ^ and (1 - e'^™'^)/{zjg'{zj)) = -1/n, then 




Now, from the Cauchy theorem, we obtain that 




1 



Finally, we conclude that 



E 



n 



1 







n . 

= — smi 




1 



f{r)dT. 



7 Conclusion 

For an arbitrary relative equilibrium of primaries in the plane, we have proved 

that each relative cquihbrium of the satclhtc generates several global branches 
of periodic solutions: for a saddle point one gets a global branch of planar 
solutions and a global branch of eight-solutions, which are truly spatial if a 
non-resonance condition holds. For a minimum point of the potential, one gets 
either two global branches of planar solutions (long and short period) and a 
global branch of eight-solutions, or only the branch of eight-solutions which is 
then truly spatial. 

A global branch may be non-admissible if the period or the norm of the 
solutions on the branch go to infinity or the branch goes to collision with one of 
the primaries. On the other hand, if the branch is admissible, then the sum of 
the jumps of the Morse indices at the critical points on the branch must be zero. 
In particular, a saddle point has to be connected with a short period minimum, 
the number of points on non-admissible planar branches is at least the number 
of saddle points and the number of saddle points on these non-admissible planar 
branches is at least one less than the number of primaries. Also, the number of 
saddle points on non-admissible branches of eight-solutions is at least one less 
than the number of primaries. 

We have applied this general result in order to describe a rather complete 
picture of the restricted three-body problem and of the restricted Maxwell ring. 

The topological degree approach, combined with the use of the orthogonality 
(or first integrals) and a systematic use of representation theory, gives informa- 
tion which is a good complement to classical analytical local calculus and allows 
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flexible applications. In particular, one may extend easily these results to dif- 
ferent potentials and to systems with more bodies. 

For concrete situations, there are many local techniques, such as normal 
form theory which often requires to check some generic assumptions ( this is 
not always done in practice), Poincare mappings, stable and unstable manifold 
decomposition of the phase space and so on. For a low dimensional bifurca- 
tion equation, there is a common starting point for these analytical methods 
and for the computation of a topological degree, that is the linearization of the 
equations. Higher order approximations may give a better local picture of the 
bifurcated solutions. But, as soon as there are resonances or more couplings, the 
analytical methods become more difficult to apply, while the topological degree 
approach can still give a complementary information on the set of bifurcating 
solutions, in particular on the global properties of the branches. It is impor- 
tant to point out that, in many relevant applied problems, one may carry out 
symbolic manipulations of high order which may be even converted into a valid 
mathematical proof using interval arithmetics. We are fully familiar with higher 
order symbolic manipulations of formal power series and the use of computer 
assisted proofs. 

With these considerations in mind, we have several papers in preparation on 
bifurcation of the whole arrangement of primaries, either as relative equilibria 
or as periodic solutions. For instance, in the case of the Maxwell ring, one 
gets n global branches of periodic solutions, each with different symmetries 
and where the central mass plays an important role, for the existence of these 
periodic solutions. Similar results were obtained for vortices, filaments, charged 
particles and nonlinear oscillators. See [9]. 
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